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Abstract 

The use of gas expansions to generate atomic or molecular beams has become a 
standard technique in nuclear and hadron physics for the production of polarized 
ion beams and gas targets. A direct simulation Monte Carlo method was used to 
understand the processes occurring in an expansion of highly dissociated hydrogen 
or deuterium gas at low densities. The results were verified in several measurements 
including time-of-flight and beam-profile determinations which showed that the su- 
personic gas expansions can properly be described by the Monte Carlo calculations. 
Additionally a new method of beam formation, the hollow carrier jet, was tested 
under the conditions of the atomic beam source operation. 
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1 Introduction 



Polarized atomic beam sources (ABS) as described, e.g., in Refs. [1] and [2] 
are utilized to provide nuclear-polarized atomic hydrogen (H) and deuterium 
(D) beams. Molecular H 2 or D 2 gas is dissociated and the essentially atomic 
gas then expands through a cooled nozzle into the vacuum. A beam of high 
brightness is then formed by a skimmer and a collimator of dimensions and 
positions to adapt the beam to the relatively small acceptance of the sub- 
sequent system of sextupole magnets. Based on the Stern-Gerlach principle, 
these magnets focus (defocus) atoms with electron-spin projection +1/2 (-1/2) 
along the magnetic field within the magnet bores. The electron-spin polarized 
beam then enters an rf transition unit, which allows to change the nuclear po- 
larization by inducing transitions between the hyperfine states. Descriptions 
of the HERMES polarized ABS are given in Refs. [3] and [4] . Details about its 
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13.56 MHz rf dissociator are found in Ref. [5], whereas the 2.45 GHz microwave 
dissociator and the atomic beam test stand are described in Ref. [6]. 

To achieve a high output intensity, the atomic beam generated by the expan- 
sion has to fulfill several requirements like high flow rate, low transversal and 
longitudinal temperature, and a high degree of dissociation. The latter can 
only be achieved, if the recombination is low. This request sets a limit to the 
gas pressure in the dissociator volume of a few mbar. The pressure within 
and at the exit of the nozzle corresponds to the transition region between 
laminar and molecular flow. There, the use of continuum-flow models (e.g., 
those based on the Navier-Stokes equations) is of restricted validity. Thus, a 
direct simulation Monte-Carlo (DSMC) method [7] was used to describe the 
processes during gas expansion. Time-of-flight (TOF) and, with the use of a 
novel type of monitor [8] , beam-profile measurements were performed to check 
the validity of the results of the simulations. The investigation of the thermal 
properties of the gas in the nozzle region gave new insights. The achieved re- 
sults could be used to improve the beam formation and to properly determine 
the beam parameters. The use of an over-expanded carrier jet, surrounding 
the inner atomic beam, had been proposed [9] to further increase the atomic 
beam intensity. This method of beam formation has been studied experimen- 
tally for the first time. The results could be interpreted by DSMC calculations. 
A detailed description of the performed work is found in Ref. [10]. 



2 Supersonic Gas Expansion 



A free-jet atomic or molecular beam can be produced by a supersonic gas 
expansion from a high-pressure gas source into a low pressure background. 
Fig. 1 shows the structure of a free expansion under continuum (steady state) 
conditions. The source is a short conical nozzle. As a result of the pressure 
difference po — Pb, the gas is accelerated. At the nozzle exit the flow may reach 
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Fig. 1. Continuum free-jet expansion of a gas into a region of background pressure 
Pb starting from a negligibly small velocity at the stagnation state, described by po 
and Tq (figure taken from Ref. [11]). 
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sonic speed, if the ratio Po/ph exceeds the critical value [11] 



G =( 1 f) ' « 

From the possible values of 7, defined as 7 = c p /cy with the heat capacities 
c p and c v , it follows that for all gases G is lower than 2.05. For p /pb < G the 
flow exits subsonically. A supersonic beam has two important characteristic 
properties. First, the velocity of the beam v increases during the expansion. 
Second, the beam parameters in the zone of silence (Fig. 1) are independent of 
boundary conditions (walls, p^), which is caused by the fact that information 
propagates at the speed of sound, whereas the gas moves faster. Therefore 
a skimmer is placed inside this zone to extract a supersonic beam. If the 
background pressure pb is small enough, a smooth transition to molecular 
flow occurs and no shock structures emerge. The beam is only affected by 
residual gas scattering. 



2.1 Thermodynamic analysis 



Considering an ideal expanding gas without viscous and heat-conduction ef- 
fects, the energy equation [11] 



v 2 



h + — = h (2) 

holds, where h is the total or stagnation enthalpy per unit mass and v is the 
mean velocity in beam direction. For ideal gases (dh = c p dT) and constant 
heat capacity c p = (7/(7 — 1))(A;b/^) one gets the maximum or terminal 
velocity (for T <T after the expansion) 



with the particle mass m and the Boltzmann constant /cb- For H2 (7 = 5/3) 
at T Q = 100 K, t>oo is 1436 m/s. For mixtures of ideal gases, an average heat 
capacity 



^- h ± 



can be used, where Xi is the fraction of the respective species. In the continuum 
limit the mean velocities of the species tend to be the same [11]. 



For isentropic expansion of an ideal gas, eq. 2 allows to deduce 

T 'l + I^M" (5) 
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Fig. 2. Free-jet on-axis properties versus distance from the source (given in source 
diameters) for a monoatomic gas (7=5/3). The temperature T, the density n, and 
the binary hard-sphere collision frequency v are normalized by the source-stagnation 
values To, no, and vq (figure taken from Ref. [11]). 

with the assumption of constant c p . Here M = v/c is the Mach number with 

the speed of sound c = ^k B T /m. Furthermore, the mean beam velocity of 
the expanding gas as function of the Mach number is derived as 



m 



7-1 2 N -V2 

1 + - M 



(6) 



With the assumption M — 1, which is expected around the nozzle exit, for H 2 
v is 718 m/s. In general, with the use of M as calculated [11] for axisymmetric 
expansion, all thermodynamic variables for the free-jet expansion can be given 
as shown in Fig. 2 for their dependences along the centerline of the expansion. 

Energy and enthalpy considerations [12], based on Eqn. 2, for an expanding 
beam of a monoatomic gas (c p = (5/2)(/c b /'/ti)), lead to the balance equation 
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-k B T = -mvl + -k B T + k B T. 



(7) 



In the term on the left side T is the stagnation temperature of the gas in the 
source, i.e., for thermal equilibrium the nozzle temperature T nozzle . The sum of 
the first and second term on the right side gives the total beam energy ^beam 
after transition to molecular flow. The third term is the energy k B T = pV, 
stored in the gas and forcing it to expand. For T = T nozzle one expects 
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2.2 Statistical definitions for the Monte Carlo method 



The Monte Carlo method is another way to obtain the parameters of the 
atomic or molecular beam by simulating binary hard-sphere collisions between 
the particles of the expanding gas. The thermal velocity v' of a particle of 
velocity v p in a beam of mean velocity v is defined by 

v p = v P - v. (9) 
The scalar pressure is defined as [7] 

P=^nm(v^), (10) 

where n is the particle density, m the particle mass, and v' p = |v p |. With 
the ideal gas law p = nk B T ', the average kinetic energy associated with the 
translational motion becomes 

E tI = ^m^=h B T tI . (11) 

For every component j this can be written as 

mvgj = k B T tIjj , j = x,y,z (12) 

For molecules the rotational and vibrational excitations have to be considered. 
They can be ascribed to an internal energy 

E int = ^CfeTint, (13) 

where ( is the number of internal degrees of freedom. 

For a gas in (local) thermal equilibrium, the fraction of the particles which 
are found within a velocity space element dv p is given by the thermal velocity 
distribution [7] 

= ( JL)'eM-l3 2 (v?, x + < ¥ + <J]d< K d< y d< a , (14) 



where (3 = \J m / (2k B T tr ) . The fraction of particles with a velocity compo- 
nent in j direction within the velocity range v' p j... v' j+dv' j, irrespectively 
of the magnitude of the other components, is obtained by integrating over 
the two other components. The distribution function for the thermal velocity 
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component then is 



exp[-/3?(u, 




(15) 



where /3j = y m/(2/c B Tt r j). The most probable thermal velocity of each com- 
ponent is zero. 



3 Simulation Program and Experimental Setup 

3. 1 The direct simulation Monte Carlo program 

The direct simulation Monte Carlo method (DSMC) [7] is a technique for 
the computer modeling of a real gas by some thousands or millions of simu- 
lated particle trajectories. The velocity components and position coordinates 
of these particles evolve in time as the particles are concurrently followed 
through representative collisions and boundary interactions in the physical 
space. The decoupling of the motion and collisions of the particles over small 
time steps and the division of the flow field into small cells are the key com- 
putational assumptions associated with the DSMC method. The time step 
should be much smaller than the mean collision time and a typical cell dimen- 
sion should be much smaller than the local mean free path. The program [13] 
has a flexible system for the specification of the flow geometry. For our purpose 
the geometry of the beam forming elements nozzle, skimmer, and collimator 
are implemented as boundary walls with temperature T for an axially sym- 
metric flow. Fig. 3 shows these elements together with the regions into which 
the flow field is divided. These regions are divided into the small cells men- 
tioned above. Additionally, the applied input flows are indicated. The main 
flow is the hydrogen flow through the cooled nozzle. In addition flow losses 
have to be compensated. So-called Specified Flows, have to be included to 
simulate the chamber pressure since particles disappear which pass the outer 
boundaries of the regions 3 and 6. The parameters of the gas, to be calculated 
for every cell of these regions, are collected in the Appendix. 

3.2 Experimental setup and data analysis 

An atomic beam test stand (ABT) had been set up [6] and equipped with 
several diagnostic devices (Fig. 4). It consists of a four-stage differentially 
pumped vacuum system with a nominal pumping speed of 14000 //s. A disso- 
ciator mounted on the first chamber produces atomic hydrogen or deuterium. 
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Fig. 3. The geometry employed in the simulation of an expansion in the ABS. 
The nozzle, skimmer, and collimator were implemented with the respective tem- 
peratures, the space in-between is divided into regions. Due to symmetry, only the 
upper half is calculated in the simulations. Specified Flows compensate for losses of 
particles at the outer boundaries of regions 3 and 6. There is no Specified Flow in 
region 7 to allow extraction of the beam parameters. 
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Fig. 4. The atomic beam test stand with the pumping stages I-IV. The pumping 
speeds are the nominal speeds for hydrogen. The signals from the quadrupole mass 
spectrometer are read out either with a secondary-electron multiplier or a Faraday 
cup. 



Two types of dissociators were used: (i) a radio frequency dissociator (RFD) 
[5], consisting of a LC circuit as a field applicator and a water cooled Pyrex 
discharge tube, and (ii) a microwave dissociator (MWD) [6], based on a plasma 
source which couples a 2.45 GHz surface wave to the discharge in an air-cooled 
Pyrex glass tube. The H (D) gas expands through a nozzle, cooled by a 120 W 
closed-cycle helium cold head. Together with a skimmer and a collimator a 
high-brightness beam is formed. This beam is analyzed by means of several 
devices. A calibrated compression tube is used to perform absolute measure- 
ments of the beam intensity. The degree of dissociation could be determined 
with the quadrupole mass spectrometer (QMS). A chopper, combined with a 
lock-in amplifier, was used to separate the beam signal from the residual gas 
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background. The intensity- related degree of dissociation a, defined as 



a = 




(16) 



was determined from the QMS signals S a and S m for atoms and molecules, 
respectively. In this equation S* = S a — 5 dlss S m is the atomic signal corrected 
for dissociative ionization in the QMS with <5 dlss = 0.0125 (0.0085) for hydrogen 
(deuterium) [6]. Furthermore, « ion = <°7<° n = 0.64 ±0.04 [14] is the ratio of 
the ionization cross sections and K det = w^/w^ 1 is the ratio of the detection 
probabilities of the QMS either with the secondary-electron multiplier (SEM) 
or with the Faraday cup (FC). Since one can assume w^ c = and therefore 
kfc = 1, /«sem can be determined from the ratio 



by measuring alternately the atomic and molecular QMS signals with the 
Faraday cup and the secondary-electron multiplier. With an accuracy of about 
1 %, «sem was determined as 0.78 ... 0.91 depending on the voltage, applied 
to the SEM. Finally, k v = v a /v m regards for the different dwell times of the 
atoms and molecules in the ionization volume of the QMS. 



3.2.1 Velocity analysis 

The velocity distribution of the particles in the atomic beam is determined 
with the time-of-flight (TOF) method. A fast chopper cuts a small bunch of 
particles out of the beam, and their arrival time at the QMS is measured. Since 
the particles have different velocities, a TOF distribution F(t) is measured. To 
keep the influence of the opening function of the chopper small, the rotational 
frequency of the motor had to be high and stable. For this reason a special 
stepping motor was used, installed in chamber III of the setup (Fig. 4). A light- 
barrier signal defined the zero point of the TOF distribution and triggered the 
oscilloscope used to store the QMS signals S(t). 

The relation between the measured TOF distribution F(t) and the velocity 
distribution f(v) of Eq. 15 with v = l cq /t is 



where l cq is the distance between chopper and QMS. The measured signal 
distribution, corrected for an offset Pq, results from the convolution of the 
TOF function F(t), the opening function of the chopper wheel A(t), and the 
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(17) 




(18) 
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response function of the electronics K(t) as 



S(t) -p* = f K{t - A) [ X F(X - ^)A(^)d^d\. (19) 
Jo Jo 

The opening function A(t) was measured at a low rotational frequency of the 
chopper wheel. At 300 Hz, used in the experiments, for the slit width of 2 mm 
and the tangential slit velocity of 75m/s, a good approximation is the half 
period of A(t) = A sm(nt/T) with A = 0.278 and T = 0.055 ± 0.004ms. 
In the response function of the electronics, K(t) = (l/r e ) exp (— t/r c ), r c = 
0.21 ± 0.01ms is calculated from the input resistance of the oscilloscope and 
the cable capacity. 



The measured signal distributions were fitted by functions 

\ 2" 



BUS * , Pi I m Mcq 



P3 



(20) 



The functions S(t) — Pq then were de-convoluted, following the method of 
Ref. [15], by calculating 



^n) = ^E(-ir 



i=0 



S(X) + r, 
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[21) 



J X=t n -iT 



for an appropriate sequence of times t n , covering the time range of the TOF 
measurement. The obtained distributions F(t n ) then were fitted by the TOF 
functions 



F(t) 



Pi 
t 2 



cxp 




(22) 



The fits yield pi as a scaling factor, p 2 according to Eqn. 15 yields the transla- 
tional beam temperature T tr x , and p^ is equal to the mean velocity v x in beam 
direction, i.e., the beam parameters to be determined. Fig. 5 as an example 
shows a measured signal distribution and the derived TOF function. 



3.2.2 Beam-profile monitor 

A beam-profile monitor [8] was used to measure intensity profiles of the atomic 
hydrogen or deuterium beam. A frame, carrying 32 gold-plated tungsten wires 
of diameter d w — 5 /xm (16 wires each in horizontal and vertical direction) was 
placed between nozzle and skimmer. At a wire spacing of 2 mm, the wires 
covered the full cross section of the beam. Beam atoms, which hit a wire, get 
sticked at the gold cover with a probability w^ u = 0.5 ±0.1 [16], i.e., they are 
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Fig. 5. Example of the measured signal distributions with the offset p$ and the 
fitted function S(t) — Pq of Eqn. 20 (left-hand side) and the TOF distribution F(t) 
(right-hand side), resulting from deconvolution, with the fitted function of Eqn. 22. 

not immediately re-emitted. If sticked, the atoms recombine with atoms of the 
surface-covering atomic layer to form molecules, which then are emitted due to 
their much lower binding energy. The recombination probability w r of sticked 
atoms thus can be assumed as 1.0. The recombination energy E r = 4.476 eV 
per H 2 molecule [17] leads to a local heating and a differential increase dR of 
the wire resistance. For a wire of length L extending in z direction, dR/dz 
is a function f(dP/dz) of the differential energy deposition dP/dz. The total 
wire resistance is 



R = R + j (dR/dz)dz = R + j f(dP/dz)dz, 



(23) 



where Ro is the wire resistance without beam. It can be assumed that (i) the 
wires are homogeneous along their length and that (ii) the locally produced 
heat is spread by radiation cooling only, i.e., that the heat transfer along the 
wire can be neglected. Under these assumptions, the function / can be deter- 
mined for each wire by application of a series of currents / and measurement 
of the voltages U. Here one has P = U ■ I, dP/dz = const = P/L, and 



R = R + J f(P/L)dz = R + L- f(P/L). 



(24) 



The measured response functions Ri(P) for the 32 wires are shown in Fig. 6. 
With some scattering, all the curves show the same behavior, an approxi- 
mately linear increase from R i at P = and a saturation towards higher P 
values caused by the increase of radiation cooling due to the increasing wire 
temperature. 

The monitor wires are positioned at different perpendicular distances yi from 
the atomic beam axis. The changes of the total wire resistances depend on the 
total flux of atoms, hitting the wires. Thus, on the one hand the measurement 
of the set of Ri — R^q yields data on the intensity distribution of the beam and 
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Fig. 6. The response functions of the 32 wires of the beam-profile monitor, measured 
by loading the wires with sequences of currents and measuring the voltages. 

its development with the distance x from the nozzle exit. On the other hand, 
the measured data can be compared with those, which result from the DSMC 
calculations. These yield atomic density and velocity distributions, n(x, y, z) 
and v(x,y,z), which can be used to calculate the differential recombination- 
energy deposition in the wires as 

dPi(x, y, z) = d w - u> s Au • w r ■ E T ■ -n(x, y, z) ■ v(x, y, z)dz. (25) 

With the use of Eq. 25 and the measured response function fi(dP/dz) = 
fi(P/L) by summation over the wire length the increase of the wire resistance 
can be calculated and compared with the measured values to check the validity 
of the calculations. 



4 Simulation Results and Experimental Verification 

4-1 Molecular hydrogen beam from a cooled nozzle 

First of all, a Monte Carlo simulation of an expansion of molecular hydrogen 
was performed. The geometry was chosen as shown in Fig. 3, the dimensions of 
the setup are given in Fig. 7. The parameters for the stream input were the gas 
temperature T gas = 300 K and the particle density no which was determined 
from the measured pressure po in front of the nozzle. The calculated flow rates 
through the nozzle, Q, were in good agreement with the measured primary 
gas flow. The parameters for the Specified Flows at the outer boundaries 
(Fig. 3) were obtained from the residual gas pressures pb measured at the 
atomic beam-test stand at the respective flow rate Q. The calculated particle- 
density distribution of the entire simulated space is shown in Fig. 7. The 
density near the nozzle follows a cos(0) distribution. Skimmer and collimator 
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Fig. 7. The calculated particle-density distribution of a molecular hydrogen beam, 
expanding in the system of nozzle, skimmer, and collimator (diameters: nozzle throat 
2.0 mm, skimmer top 6.4 mm, collimator 11mm; distances: nozzle to skimmer-top 
15 mm, skimmer-top to collimator 50 mm; T QOZZ \ c 100 K; primary molecular gas flow 
1 mbar^/s). 

form a low-diverging and sharp-bound molecular beam. The density at the 
beam edge decreases by two orders of magnitude over a distance of 1 mm at a 
beam diameter of about 15 mm. 



Fig. 8 shows the calculated on-axis beam properties. Within the nozzle (x = 
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Fig. 8. On-axis beam properties, calculated for a H2 beam in the setup of Fig. 7 
as function of the distance to the nozzle exit (n: particle density; Tt r : translational 
temperatures describing the velocity distributions in x (full circles) , y and z direction 
(open circles); Ti nt : rotational temperature; v x : mean particle velocity in x direction 
(along the beam axis)). The pronounced variations of T tr , Tj n t, and v x between 
x = 2 mm and x = 65 mm are caused by the sampling of the DSMC computer code 
as discussed in the text. 
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— 15 mm to mm), T tr and T- mt drop to the temperature of the nozzle, i.e., the 
gas is thermalized, the particle density n rises, and the mean velocity v x is 
low. At the exit of the nozzle (or shortly before) the expansion into the low- 
pressure region starts. The particle density decreases with increasing distance 
from the nozzle. 



The increase and drop of T tr and T int as well as the drop and rise of v x between 
x = 2 mm and x = 65 mm, the position of the collimator, are due to the sam- 
pling by the DSMC computer code. Every data point of the distributions of 
Fig. 8 results from averaging over the molecules in the entire cell. It includes 
beam molecules as well as residual gas molecules that will not pass the col- 
limator. Genuine beam parameters can only be extracted from the cells with 
x > 65 mm, since no Specified Flow is applied to region 7 (Fig. 3). Regarding 
this fact, one can state that for x > 2 mm T tr , T; nt , and v x essentially stay 
unaffected. The reason is that the collision rate is too low to enable thermal 
equilibrium between translational and internal energy, which requests about 
100 bounces per molecule. The internal (rotational) temperature "gets frozen" 
in the region of the nozzle exit. Contrary to Tj nt , T tr shows a pronounced de- 
crease between x = 0, the nozzle exit, and x ~ 2 mm, the position of the 
"freezing surface" . A few bounces only are required to thermalize the transla- 
tional degrees of freedom. In the same x interval, the mean velocity in beam 
direction, v x , increases with the decreasing translational beam temperature. 
Thermal energy is converted into directional motion. The temperatures T tr in 
y and z directions are lower than that in x direction, because particles with 
higher transversal velocity leave the beam center, and thus the temperature 
on the axis drops. This effect may be denoted as "geometrical cooling". 



In the upper part of Tab. 1, the values for the mean velocity v x and the 
translation temperature T trjX are compared as they result from measurement 
and DSMC calculation, both performed under the boundary conditions given 
in Fig. 7. The calculated mean velocity is close to the measured one, but the 
resulting temperature is appreciably too high. The problem was studied [18,19] 
and it was found that the temperature contrary to the mean velocity, strongly 
depends on the parameters of the collision processes used in the simulation 
code. It seems that certain parameters, chosen by the editors of the simulation 
code [13], are not fully correct for the low-temperature region. Furthermore, 
the influence of background gas in the nozzle-skimmer region over the beam 
properties was investigated [18,19]. It could be shown that by the modification 
of parameters, essentially of scattering cross sections, the beam temperature 
can be adopted to the measured one, while the mean velocity stays almost 
unaffected (Tab. 1). 
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Table 1 

Comparison of the calculated and measured parameters of (i) a pure molecular and 
(ii) a partly dissociated hydrogen beam. The errors of the measurements are statisti- 
cal and systematical, resulting from those in the determination of T and r e (section 
3.2.1), respectively. The asymmetric systematical errors of the simulation results 
originate from discrepancies of the simulation distribution and the Maxwellian dis- 
tribution used for the analysis of the experimental data [20]. The boundary condi- 
tions are given in Fig. 7. The measured mean velocity in beam direction, v x , and the 
temperature T tr)X result from fits according to Eqn. 22. M x = ^ x /\/7^B^tr,x/ m is 
the axial Mach number and Voo is the maximum beam velocity according to Eqn. 3. 
The simulation calculation of Refs. [18,19] were initiated by the discrepancy en- 



countered with the original parameters. The dej 
partly dissociated beam was a = 0.63. 


;ree of dissociation (Eqn. 16) of the 




v x (m/s) 


T tr , x (K) 


M x 


Voo (m/s) 


molecular beam 










measurement 


1274 ± 8 ± 13 


19.0 ± 1.1 ±0.9 


3.52 




simulation (original) 


1334 ± 12 ±[ 5 


33.3 ± 1.6 tl 


2.79 


1436 


simulation (Refs. [18,19]) 


1371 ± 2 ± 5 15 


19.0 ±0.2 ±° 3 


3.79 


1436 


partly dissociated beam 










atoms 










measurement 


1750 ± 47 ± 24 


25.7 ± 4.9 ± 1.5 


2.94 




simulation 


1760 ± 20 ±f 9 


41.0 ± 3.4 ±° 3 


2.34 


2031 


molecules 










measurement 


1579 ±51 ± 17 


23.7 ± 7.1 ± 1.0 


3.91 




simulation 


1590 ± 33 t\ 7 


44.0 ± 4.3 ±1 


2.89 


1436 



4-2 Partly dissociated beam from a nozzle 

With the use of the QMS, the TOF distributions of the atoms and molecules 
in a partly dissociated hydrogen beam were measured. The measurement was 
performed with the same setup of nozzle, skimmer, and collimator which was 
used with the pure molecular hydrogen beam. The dimensions were those given 
in Fig. 7. The degree of dissociation a, entering the following calculations, 
was between 0.50 and 0.80 depending on the experimental conditions. TOF 
distributions were measured with nozzle temperatures T QOZZ \ c = 70, 100, 150, 
and 200 K and primary molecular flows Q between 0.5 and 6.9mbar//s. The 
mean velocity and beam temperature for the atomic and the molecular fraction 
of the beam, measured with T QOZZ \ c = 100 K, Q = lmbar//s and a = 0.63 and 
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resulting from the fits according to Eqn. 22, together with the Mach numbers 
are given in the lower part of Tab. 1. 

The temperature of the gas in the plasma of the microwave dissociator [6], 
used in these measurements, may reach values of T plasma = 2000 K [21] or 
even higher. Since the plasma end is near to the nozzle entrance, one has to 
investigate, whether the gas reaches thermal equilibrium in the nozzle before 
the expansion. According to Eqn. 8, the beam energy (in units of /cB^nozzie) 
is expected to lie in the range 3/2 to 5/2 for full temperature equilibrium of 
the gas in the nozzle. The left-hand part of Fig. 9 shows four sets of beam 
energies, calculated from v x and T tr x of both species H and H 2 , as function of 
the primary molecular flow Q. 

For Q < lmbar//s, the beam energies lie in the expected range 3/2 to 5/2. 
For higher input flows, however, it is higher than expected. The discrepancy 
increases with decreasing nozzle temperature T nozz i c and increasing gas flow Q. 
This finding would be explained by the assumption that the temperature of 
the gas in the nozzle, T , is higher than T nozz i c , i.e., that the gas does not reach 
the thermal equilibrium. For Q = 1 mbarZ/s and T nozzle , this explanation was 
studied by DSMC calculations. The right-hand part of Fig. 9 shows the calcu- 
lated on-axis temperature of the gas in the nozzle cone (T nozz i e = 100 K) as a 
function of distance to the nozzle exit for four temperatures Tpi^^ of the gas 
entering the nozzle. For Tpi asma = 300 K, the thermalization is complete, while 
for higher T Plasma the thermalization is incomplete, T > T nozz \ e , and hence the 
beam energy exceeds 5/2 /c B T nozz i c . Because the exact gas temperature in the 
plasma is not known, in the further DSMC calculations T p i asma = 3000 K was 
used as the temperature of the gas, entering the nozzle. 
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Fig. 9. Left-hand side: beam energies, as defined in Eqn. 7 and calculated from the 
measured velocity distributions of H and H2 , as function of the primary molecular 
gas flow Q. Each of the four sets is given in units of kBT nozz \ e with the respec- 
tive nozzle temperature. Right-hand side: on-axis gas-temperature distributions in 
the nozzle (T nozz i e = 100 K), calculated with the DSMC code for Q = ImbarZ/s, 
a = 0.67 and four temperatures of the plasma in the microwave dissociator, as 
function of the distance to the nozzle exit. 
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Fig. 10. The calculated density distribution of the atomic (upper plot) and molec- 
ular (lower plot) fraction of a partly dissociated hydrogen beam (a = 0.63) for 
^nozzle = 100 K and Q = ImbarZ/s. 



The calculated density distribution of the atomic and molecular fraction are 
shown in Fig. 10. As for the pure molecular beam, low-diverging atomic and 
molecular beams with sharp boundaries are obtained, formed by skimmer and 
collimator. The distributions of the on-axis beam properties (particle density, 
translational and internal temperature, and mean velocity) are similar to those 
shown in Fig. 8 for the pure molecular beam. The calculated mean velocities 
and temperatures for the atomic and molecular fraction in Tab. 1 are juxta- 
posed with the measured values. Both the measured and calculated values of 
the mean molecular velocity are in good agreement. They show that in the ex- 
pansion of a gas mixture the heavier species (molecules) are accelerated above 
the theoretical limit of a pure expansion, which would yield t> Xj H 2 = ^x.h/v^- 
As for the pure molecular beam, one finds a deviation of the simulated from 
the measured temperatures, while the mean velocities agree. The reason for 
this difference is the same as discussed for the DSMC calculations for the 
molecular beam. 
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Fig. 11. Calculated distributions of the atomic hydrogen density n(x, y, z = 0) (left- 
-hand histogram) and of the mean velocity of the atoms v(x,y,z = 0) (right-hand 
histogram). The data of the atomic beam stem from a DSMC calculation of the 
expansion of a partly dissociated hydrogen beam (rf dissociator with plasma tem- 
perature T p i asma = 3000 K, degree of dissociation a = 0.76), and nozzle temperature 
^nozzle = 100 K. Because of cylindrical symmetry the dependence on z at y = is 
identical to the y dependence at z = shown. 



4-3 Beam profiles 



The data, obtained by the DSMC calculations for a partly dissociated hydro- 
gen beam, could be further compared with the results of beam-profile measure- 
ments, made with the beam-profile monitor described in the section 3.2.2. The 
distance between nozzle and skimmer had to be increased from 15 mm (Fig. 7) 
to 50 mm to allow positioning and moving of the monitor in-between along the 
beam axis. Therefore, only the expansion from the nozzle exit into the low- 
pressure region between nozzle and skimmer was considered. The input geome- 
try of the DSMC code was reduced to this region and the skimmer was assumed 
to be the transition to vacuum. The calculated density distribution n(x, y, z) 
and velocity distribution v(x, y, z) = ^Jv^.(x, y, z) + v^(x, y, z) + v%(x, y, z) are 
shown in Fig. 11. The left-hand histogram indicates, how the total density 
of particles decreases from the beam axis (y = 0) to zero. Contrary to that, 
the mean velocity of the atoms, shown in the right-hand histogram, decreases 
only slightly from the on-axis value of about 1700 m/s, comparable to those 
of Tab. 1, to about 1400 m/s. In the covered range of x, the density shows a 
distinct decrease, whereas the velocity stays more or less unchanged due to 
the fact that already the minimum x position of 10 mm is far behind that of 
the freezing surface at about 2 mm. 

The beam monitor allows to measure the change of the resistance of each of 
the wires by surface recombination of the hydrogen atoms in the beam (cf. 
section 3.2.2). In the used coordinate system, given in Fig. 1, each of the wires 
extends along the z direction with a perpendicular distance y to the beam axis 
in a plane, positioned perpendicular to the beam axis at distance x from the 
nozzle exit. With the calculated n(x,y,z) and v(x,y,z), the differential heat 
deposition in a wire can be calculated with the use of Eqn. 25. Taking into 
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Fig. 12. Calculated (left-hand part) and measured (right-hand part) distribution of 
the wire resistances (the geometry and the experimental parameters are those given 
in Figs. 10 and 11). 



600 ^ + 

550- ♦ frM| , 

500- t | ^ 
450 ' 

-10 10 -10 10 

y (mm) y (mm) 

Fig. 13. Comparison of the calculated (circles)and measured (triangles) distributions 
of the wire resistances at distances between wire plane and nozzle exit of 10 mm 
(left-hand side) and 20mm (right-hand side). The lower-lying distributions (full 
symbols) were obtained for the geometry and the experimental parameters given 
in Figs. 10 and 11, whereas the higher-lying distributions result with a primary 
molecular flow Q enlarged from lmbar//s to 2mbar//s. 
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account the measured response functions of Fig. 6 and summation over the 
wire length yields the distribution R(x, y). The density of the hydrogen atoms 
at the wire ends, n(y = ±25 mm), is very small compared to that of the beam 
axis. Thus recombination effects on the frame could be neglected. In Fig. 12, 
the calculated distribution of the resistances is compared with the measured 
one. Both distributions are not smooth due to the variations in the response 
functions of the wires. Good agreement, however, is found for the shape of the 
distributions and their absolute values. In more detail, this is demonstrated in 
Fig. 13. There distributions are shown for two distances between wire plane 
and nozzle exit, x = 10 mm (left-hand part) and x = 20 mm (right-hand 
part). The lower-lying distributions result for the geometry of Fig. 10 and the 
experimental parameters of Fig. 11, i.e., they are cuts through the distributions 
of Fig. 12 at x = 10 mm and 20 mm, whereas the higher-lying distributions 
are found for an increased primary molecular flow Q = 2mbar//s. The good 
agreement confirms the validity of the DSMC calculations. 
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4-4 Hollow carrier jet 

The carrier jet method was proposed [9] to increase the phase-space density of 
the atomic beam and thus to reach a higher intensity through the collimator of 
the atomic beam source. An over-expanded carrier jet, surrounding the inner 
atomic beam, was predicted to cool and to confine the inner beam (Fig. 14). 
The mixing of the two gases has to be small and the carrier gas has to be 
removed by the skimmer and pumped away. The present measurements and 
calculations were extended to investigate for the first time the idea of the 
carrier jet. A variety of inner/outer gas combinations were studied, namely 
H/H 2 , D/D 2 , D 2 /H 2 , D/He, and Ar/N 2 . The test stand (Fig. 4) had to be 
modified by replacing the turbo-molecular pumps of the first chamber by a 
roots-pump system with a nominal pumping speed of 1000 Z/s to master the 
requested [9] high carrier-gas flows (up to 40mbarZ/s). The aluminum carrier- 
jet nozzle combines a conical inner nozzle with an outer, ring-like Laval-type 
nozzle to create an outer hollow beam surrounding an inner nozzle beam 



4-4-1 Low-mass gases 

Starting with molecular deuterium D 2 as inner gas and molecular hydrogen 
H 2 as carrier gas, measurements were performed for a wide range of input 
parameters. The nozzle-throat diameter was 2 mm and the distance between 
the nozzle exit and the skimmer was 30 mm. The nozzle temperature was kept 
at T nozz i c = 100 K. For three values of the inner D 2 flow of 1, 4, and 7mbar//s, 
the outer H 2 flow was varied from up to 30 mbarZ/s to approach the suggested 
ratio [9] of about 40. The QMS signals, which give the particle density in the 
ionization volume, were multiplied by the mean velocity, determined by TOF 
measurements, to yield the on-axis intensities. The results are presented in 
Fig. 15. As the left-hand part of the figure shows, no increase of the on- 
axis D 2 beam intensity is found. On the contrary, especially the intensity 
shows a strong decrease at ImbarZ/s primary D 2 flow. The results of Fig. 15 
were confirmed in further measurements with different dimensions of the outer 
nozzle, with different nozzle-to-skimmer distances, and with H/H 2 , D/D 2 , and 



(Fig. 14). 



carrier gas 




Fig. 14. The nozzle and the principle of the hollow carrier jet. 
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Fig. 15. Measured on-axis intensities of the inner gas D2 (left-hand side) and the 
carrier gas H2 (right-hand side) as a function of the primary carrier gas flow Qn 2 
for primary inner gas flows Qd 2 = 1 mbar/ / s (triangles down) , 4 mbar/ /s (triangles 
up), and 7mbar//s (dots). The intensities Id 2 and Ih 2 are given in arbitrary, but 
identical, units and, thus, the given values can be compared directly. 

D/He as inner/outer gases. The measured high on-axis intensities of the carrier 
gases indicate a pronounced mixing of the outer gas into the inner beam. This 
is obvious, too, from the beam parameters, mean velocity and temperature, 
measured for the (inner) D 2 beam as shown in Fig. 16. The gas of the inner 
beam gets more and more accelerated with increasing carrier-jet flow. No 
cooling effect on the inner beam was measured. 

Because of the pronounced discrepancy between the measured data and the 
prediction [9], the DSMC program was used to understand the mechanisms in- 
volved. The calculations were performed for a nozzle-throat diameter of 2 mm, 
equal to that used in the measurements, and a slightly increased distance be- 
tween nozzle exit and skimmer top of 35 mm. The obtained density distri- 
butions are shown in Fig. 17. There, the upper plot shows that of the D 2 
beam of gas flow Qd 2 — 7mbar//s, expanding without carrier gas. The cen- 
tral plot presents the distribution of the D 2 fraction (Qr> 2 — 7mbar//s from 
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Fig. 16. Mean velocity v x and beam temperature 7tr,x of the inner gas D 2 as a func- 
tion of the primary carrier gas flow Qn 2 for primary inner gas flows Qb 2 = 1 rnbarZ / s 
(triangles down), 4mbar//s (triangles up), and 7mbar//s (circles). 
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Fig. 17. Calculated distributions of the density in the expansion of an inner D2 beam 
without carrier jet (upper plot) and the densities of the inner, D2, fraction (central 
plot) and the outer, H2, fraction (lower plot) in the expansion with a carrier-jet. 
The gas flows were chosen as Qb 2 = 7mbar//s and Qh 2 = 21 mbar//s. 

the inner nozzle as in the upper plot) in the expansion with an H2 carrier jet 
(Qh 2 — 21mbar//s) from the outer ring-shaped nozzle. The lower plot shows 
the density distribution of the H 2 carrier molecules. No significant difference 
is found between the distributions in the upper and central plot. In agree- 
ment with the measurements, the comparison of the central and the lower 
plot shows that the on-axis density of the H 2 carrier gas even can exceed that 
of the D2 beam. A large amount of the carrier gas passes through the skimmer 
in contrast to the predictions. It seems that the mixing coefficients of hydro- 
gen/deuterium are too high at these densities to create the carrier jet effect. 
The Monte Carlo simulations again verified the experimental results. It can 
be concluded that the Navier-Stokes equations, applied for the predictions [9] 
are not valid in the investigated flow range. 



4-4-2 Heavier gases 

The mixing of two gases depends on their atomic or molecular diameters 
and masses. Thus an additional measurement was made using argon as inner 
and nitrogen as carrier gas. The dependences of the argon and nitrogen on- 
axis intensities on the nitrogen carrier-jet flow are found in Fig. 18. Here, 
contrary to the measurements with hydrogen and deuterium, the QMS-signal 
distributions are given without multiplication by the velocity, because the 
measured mean beam velocity v x does not change in the studied range of Q^ 2 
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Fig. 18. The QMS signal of the inner argon gas (left) and the nitrogen carrier gas 
(right) as a function of the carrier gas flow Qn 2 f° r inner argon flows Qxx = 1 mbarf/s 
(triangles), 4mbar//s (squares), and 7mbar//s (dots). The distance between skim- 
mer-top and collimator was 35 mm and T nozz i e = 100 K. 

(see Fig. 19). The distributions of Fig. 18 are clearly different from those, 
measured with D2 as inner and H2 as carrier gas. For the highest argon flow of 
7mbar//s, the measured on-axis argon intensity increases by a factor two from 
the argon expansion without carrier gas. Almost no nitrogen was detected in 
the QMS at argon flows higher than 3mbar//s. 

The confirmation of the carrier-jet effect is found in the results of the TOF 
measurements. Fig. 19 shows the parameters of the beam with inner (argon) 
flow Qxr = 7mbar//s as a function of the nitrogen-carrier flow Q^ 2 . A clear 
cooling effect could be seen with a minimum around Q^ 2 = 3mbar//s, exactly 
where the intensity maximum appears (upper curve Fig. 18 left side). Thus 
the N2 beam cools the Ar beam without mixing because the mean velocity 
stays constant. 



The density distribution calculated by the Monte Carlo simulation program 
of the argon expansion without (upper) and with carrier jet (middle) is shown 
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Fig. 19. Beam temperature T tr;X and mean velocity i> x of the inner argon gas as 
a function of the carrier gas flow rate Qn 2 for the inner argon flow of 7mbar£/s 
(distance between skimmer top and collimator 35 mm, T nozz \ c = 100 K, cf. Fig. 18). 
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Fig. 20. Calculated distributions of the density in the expansion of an inner argon 
beam without carrier jet (upper plot) and the densities of the inner, argon, frac- 
tion (central plot) and the outer, N2, fraction (lower plot) in the expansion with a 
carrier-jet. The gas flows were chosen as Qav = 7mbar//s and Qn 2 = 4mbar//s. 

in Fig. 20. The on-axis argon density is increased with the carrier gas flow 
rate Qn 2 = 4mbar//s compared to the carrier- free expansion. The flow rate 
through the skimmer increases, while the divergence of the beam is almost 
the same. The on-axis density of the nitrogen carrier gas is very small as well 
as the flow through the skimmer in contrast to that of H2 Fig. 17. In both 
cases, the results of the DSMC calculations are consistent with those from the 
measurements. The mixing of the two gases is restricted to the boundary layer, 
and it obviously is much smaller than for hydrogen into deuterium in that 
flow regime. Therefore higher densities in the expansion would be required to 
reach a carrier-jet effect for deuterium/hydrogen. The degree of dissociation, 
however, for flow rates Q > 7mbarZ/s is too low to achieve higher intensities 
by atomic beam sources. 



5 A further application of the DSMC method 



Apart from optimization of beam formation, the direct simulation Monte Carlo 
(DSMC) method can be applied to generate the input data for Monte Carlo 
simulations (see, e.g., Ref. [22]), where the atoms are tracked through the fields 
of the sextupole magnets in polarized atomic beam sources (ABS) as it has 
been mentioned in the introduction. Position and velocity distributions of the 
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Fig. 21. Models used to define the input distribution of the atoms entering the 
sextupole system in tracking calculations (a: molecular flow model, b: laminar flow 
model, c: model of a flow in the transition region). 

atoms on the collimator surface are needed as input data to achieve a correct 
calculation of the output intensity and polarization of the ABS. In addition 
to the velocity distributions in beam direction, determined by time-of-flight 
measurements like those of Ref. [23], a variety of models were used to create 
a sample of atoms passing the collimator [24]. 

• In the molecular flow model the connecting line between two randomly 
distributed points, one in the nozzle exit and the other in the aperture of 
the collimator, defines the direction of the atom (Fig. 21a). 

• The laminar flow model uses the apex S of the cone, defined by the open- 
ings of nozzle and collimator, as the first point and a random point in the 
aperture of the collimator as the second one (Fig. 21b). 

• The generalization is the model of a flow in the transition region (Fig. 21c). 
The first point is generated on a virtual nozzle and the second one is a 
random point in the aperture of the collimator. A molecularity parameter 
is defined as 

log(l - ^) 

K mol = f^, (26) 

M 1 " fe) 

where d VQ , d pn and d co are the distances of the virtual nozzle, the physi- 
cal nozzle and the collimator to the point S. For a given nozzle-collimator 
distance, the choice of K mo \ defines the position of the virtual nozzle. 

For the molecular flow model K nio \ = 1 and for the laminar flow model 
-Kmol — 0. Furthermore, K mo i > 1 means that the virtual nozzle lies between 
the real nozzle and the collimator. To decide, which of these models yields a 
satisfactory approximation of the real distributions of the particle directions 
and velocities, the model parameters can be compared with those resulting 
from DSMC calculations. 



For the flow models, the unknown thermal non-axial velocity distributions 
can be approximated by a convolution of a function describing the nozzle- 
collimator geometry and a velocity distribution resulting from the measured 
mean axial velocity v x and the axial beam temperature T tr x , 
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Fig. 22. Mean radial velocities of the hydrogen atoms in the collimator plane, cal- 
culated in the framework of the models of gas expansion for different distances y to 
the beam axis (left-hand side) and non-axial beam temperatures (right-hand side). 
The curves are given as function of the molecularity parameter K mo \. 

where j — y and z denotes the components perpendicular to the beam di- 
rection in the coordinate system of Fig. 1. The parameter Cj is given by the 
nozzle-collimator geometry and the choice of the position of the virtual nozzle, 
i.e., the value of K mo \. Furthermore, v'j max = Cj ■ v x (j = y, z) is the maximum 
transverse thermal velocity possible at the mean axial velocity v 
temperatures T tT j of the gas with the use of eq. 12 are 



and Cj. The 
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The left-hand part of Fig. 22 shows the mean radial velocity of particles cross- 
ing the collimator aperture at different radial distances y from the beam axis. 
The right-hand part shows the non-axial beam-temperature dependence. The 
curves are given as a function of the molecularity parameter K mo \, which enters 
into the Eqns. 27 and 28 via the choice of Cj. 

The parameters, obtained in the framework of the flow models, can be in- 
terpreted by comparison with those from DSMC calculations, presented in 
Fig. 23. By this comparison, the validity of the three models can be discussed 
concerning the distributions of the input data for the trajectory calculations. 

• The molecular-flow model (K mo \ = 1) with the assumption of particle emis- 
sion from the real nozzle allows T trjZ to be different from zero, because not 
all the trajectories lie within planes containing the beam axis. The non-axial 
temperatures T tr y and T tr z , resulting from the DSMC calculation with val- 
ues around 2K (Fig. 23), lie higher than the small values of ~ 0.1K from 
the model approximation. The DSMC calculation, however, may yield too 
high temperature values with the standard parameters in the code as it 
was discussed in section 4.1. Modification of the parameters might yield a 
lower temperature. The analysis of the density-distribution data, obtained 
with the DSMC code (Fig. 10) yields < K mol < 1, which supports the 
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Fig. 23. Beam parameters of an expanding, partly dissociated hydrogen beam ex- 
pansion calculated with the DSMC code (primary gas flow Qh 2 = 1 mbar//s, degree 
of dissociation a = 0.8, T DOZZ \ e = 100 K, T p i asma = 3000 K). In the lower plots full 
circles denote the x, open circles the y, and stars the z direction. 

applicability of this model. With emission of randomly distributed particles 
from the nozzle, it is the simplest model one can think of. It is thought to be 
applicable, when high accuracy is not requested in the tracking calculations. 
• The laminar flow model (K mo \ = 0) was frequently used, since the radial 
beam profile, measured behind the collimator [22], was considered to be 
sufficiently consistent with that predicted by the model. Particle emission 
from a point like source on the beam axis implies that the non-axial beam 
temperatures are zero. The results of the DSMC calculations (Fig. 23), 
however, show that the transversal beam temperatures are different from 



zero (even, when they are reduced from the calculated T, 



tr,y,z 



2 K to IK 



due to the reasons given above for the molecular flow). A temperature of 
IK corresponds to a root mean square velocity of {v'^j) 1 ^ 2 = 91m/s. Total 
omission of those particle trajectories, not lying in planes containing the 
beam axis, means neglecting atoms with an angular velocity component 
and modified trajectories in the sextupole system due to the centrifugal 
force. Therefore, the laminar flow model is thought not to be suited to yield 
the appropriate input data. 

For the model of a flow in the transition region, the transversal beam tem- 
peratures from the DSMC calculation, reduced to 1 K as discussed above, 
according to Fig. 22 yield K mo \ ~ 4. For this value, the increase of the 
transversal beam velocity v y from the beam axis (y = 0) to the aperture 
radius of the collimator (y = 5.5 mm) roughly match that of the DSMC 
data (Fig. 23). According to its definition, Eqn. 26, the virtual nozzle for 
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this value is positioned between the real nozzle and the collimator, which 
might be explained with particle emission from the freezing zone, i.e., the 
zone, where the laminar flow turns to the molecular flow. 

The molecular-flow model and the model of a flow in the transition region are 
thought to be suited to yield the input-data distributions for the trajectory 
calculations in the sextupole magnets of a polarized ABS. The latter model 
might be the superior one, when the position of the virtual nozzle can be fixed 
to the appropriate position. 



6 Summary and Conclusion 

It has been shown that the direct simulation Monte Carlo (DSMC) method [7,13] 
is an excellent tool to describe the processes occurring in the expansion of light 
and also heavier gases in the transition region between laminar and molecu- 
lar flow. The results of the calculations were confirmed by the measurements, 
performed at an atomic beam test stand (Fig. 4, [6]) with the use of an rf 
dissociator [5], a microwave [6] dissociator, and a novel atomic beam-profile 
monitor [8]. The origin of the discrepancies between simulated and measured 
temperatures was found, the problem could be solved [18,19] by modification 
of the appropriate input parameters of the DSMC code [13]. The predicted 
carrier-jet effect [9] could not be observed for hydrogen and deuterium at the 
operational parameters of atomic beam sources. It was, however, observed 
for an argon beam surrounded by a molecular nitrogen-carrier jet. Both ex- 
perimental findings are consistent with the results of DSMC calculations. The 
problems, occurring with the so far used start generators for particle-trajectory 
calculations through the sextupole magnets of polarized atomic beam sources 
were studied. Uncertainties in the used models remain, and the distributions 
from DSMC calculations are considered to present the most reliable start gen- 
erator. 
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Appendix 



The output parameters of the DSMC program [13] in the two output files 
of the first column are shown in the table below. The x- and y-coordinate 
refer to the center of the corresponding cell. The mean velocity of species % in 
j-direction is v* = Vj — Ai>j. The overall temperature T is calculated as [7] 

T = (3T tr + Crot^int* + (vibT^t)/ (3 + Crot + Cvib), 

where ( are the mean numbers of degrees of freedom of the species. It has to 
be mentioned that in the axially symmetric case utilized here, x is the beam 
direction, y the radial direction and z the circumferential direction, i.e., v z is 
tangential to the angular velocity. 



output file 


parameter 


description 


uszgr.txt 


\' /~\ f\ r\ y> T\ 

A CvJUKJJ 


x coordinate 




Y LUUKU 


y coordinate 


(for the 




mass density p 


entire 


TR TEMP 


kinetic temperature T tr 


beam) 


ROT TEMP 


temp, of rotations T[°£ 




VIB TEMP 


temp, of vibrations 




OV TEMP 


overall temperature T 




MACH 


Mach number M 




U 


mean velocity 




V 


> 


f X , Vy, V z 




w 

) 


of all species 


ds2gm.txt 


N DENS 


number density n 




TTX 


kinetic temp, in x Tt r , x 


(for every 


TTY 


kinetic temp, in y Tt r , y 


species) 


TTZ 


kinetic temp, in z Xt rjZ 




U DIF VEL 


mean velocity difference 




V DIF VEL 


> of a species i in x, y, z 




W DIF VEL 


A<, Av y , Avl 
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